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Abstract 

The problem of a moving semibounded magnetoactive plasmas in a 
plane gravitational wave is research on the basis of the self-consistent 
equations, obtained earlier by the Authors. 



1 Introduction 

In the series of articles ([T], [5], [3], [S]), based on exact stationary 1 solutions 
of the set of Maxwell equations and the equations of relativistic magnetohydro- 
dynamics for an unlimited homogeneous plasma, we established an essentially 
nonlinear behaviour of a highly magnetized plasma in the field of gravitational 
radiation propagating across the magnetic field. On reaching some critical con- 
ditions, the plasma begins accelerating up to subluminal velocities in the di- 
rection of the gravitational wave propagation. Simultaneously, the magnetic 
field strength greatly increases. This article is devoted to an analysis of the 
boundary-value problem for a semibounded plasma, with a goal to study the 
mechanism of the origin of a gravimagnetic shock wave. 

2 Self-consistent equations of motion of a mag- 
neto active plasma in the PGW field 

In Ref. from the condition of coincidence of dynamic velocities in the 

energy-momentum tensors (EMT) of a perfect fluid and the electromagnetic 
field, the equations of relativistic magnetohydrodynamics of plasma in the ar- 
bitrary gravitational field were obtained, and their exact solution was found 
in the background of the plane gravitational wave (PGW) metric, correspond- 
ing to an initial state of a homogeneous electroneutral plasma with a frozen-in 
homogeneous magnetic field. 

A general property of the above solution is the presence of a singularity on 
the hypersurface 

A{u) = 1 - a' [e^'^^") - l] = 0, (1) 



a' = sin^ n/ATT{so + po) (2) 

is a dimensionless parameter while /3(u) is an arbitrary function of the retarded 
time u (the PGW amplitude). This singularity in a weak PGW (/? <C 1) is 
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reached under the condition: 

2/3{u)a'^ ^ 1 . (3) 

An analysis of the solution shows that, when the condition 1^ holds, a shock 
wave is formed in the plasma, propagating with a subluminal velocity in the 
direction of PGW propagation. A necessary condition for the excitation of 
a shock magnetohydrodynamic wave (GMSW) is a high magnetization of the 
plasma: 

> 1 . (4) 

However, the above solution of the RMHD equations is essentially stationary 
(it depends only on u, the retarded time) and corresponds to an initially ho- 
mogeneous magnetoactive plasma. Therefore the solution obtained in Ref. 0]] 
cannot describe the dynamics of the shock wave excitation mechanism. 

In this paper, we will study the problem of PGW distribution in an isotropic 
magnetoactive plasma with boundary conditions, supposing that the null hyper- 
surface: 

So : It = (5) 
is the surface of PGW front, i.e., the PGW is absent at u ^ 0, - 

/3(u)|„^o = 0; /?'(u)|„^o = 0; L(u)|„^o = l- (6) 

In the absence of a PGW, the plasma is at rest, i.e., 

'0(u,w)|„^o = "00(2;^) ; (7) 

P|«^o =Po(a;^) ; £\u!^o ^ £o{x'^) , (8) 

where ipo ,pq ,eo - are some given functions of the variable x^; 

Thus, due to ©: 

- {dutp)\ui:o = (3«0)|«s:o • (10) 

In the absence of a PGW, the energy density, pressure and magnetic field 
strength are identical everywhere: 

e|usSo = £o = Const ; p|„^o = Pa = Const ; (11) 

H\u^a = Ho = Const ^ i'^^^o = -^Ho{v - u), (12) 
at w ^ ; V w en. 

In Ref. [4] it has been shown that the plasma macroscopic parameters are 
the following: 

duZ I a„z 
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ij2 = -L-'^e'^'^duZd^Z , 



where 



£ = L-hoy/-duZd„Z 
V2^ 



Z 



Ho 



(14) 
(15) 

(16) 



is a dimensionless function. 

The function Z{u,v) should satisfy a boundary condition on the spacehke 
hypersurface J^r ■ u — v {x — 0).) Let us put, following Ref. [1], on this 
hypersurface: 

H\u=v = Const. (17) 
Then, taking into account the relation (see Ref. [I]): 



we bring the boundary condition (jl7p to the form: 



iduZd,Z\^^^ 



1, 



(18) 



(19) 



Taking into account these conditions, the set of RMHD equations is reduced 
to a single differential equation: 

a..Z + /3UZ=^$(s+p)(^')' 



de' 



^0 



dvvZ 



dv.ii Z 
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duZ 
duvZ 



{d^zf [duzy duZd^z 



dp 

de 



dyyZ 



duuZ , 2 duvZ 



_{dyZf (a„z)2 duZdyZ 

which for a non-relativistic plasma becomes 



(20) 



duvZ + P'd,Z 



dvv Z 



-^e-^^^-d.,ZduZx 

, (^uuZ 2 ^uvZ 



{dyZf {duZf ^duZdyZ 
(for details see Ref. [S]). 



(21) 



3 Search for a solution in a weak gravitational 
wave 

Since the GMSW is formed in a highly magnetized plasma, in what follows we 
shall suppose that a ([2]) is a large parameter (see the condition Q). Therefore, 
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Eq. (|2ip . being expanded in the small dimensionless parameter of the problem 
a^^, takes the form: 

duvZ + /3'dyZ ^0 (22) 

Solving it, we obtain the first approximation, satisfying the conditions indicated 
above: 

Zo = {v~u)e-''^ -2 smh{l3{u))du (23) 

"'0 

The magnetic field behaviour corresponding to the solution is shown in 

Fignj. 




Figure 1: Dependence of the magnetic field inten- 
sity on the retarded time calculated according to the 
solution (f23)) . 

To verify the linear approximation applicability, we substitute the solution 
found to the r.h.s. of Eq. ([2T|) and obtain: 

duvZ + l3'd,Z = ^u,v), (24) 

where $(m, u) is a very unwieldy expression. It is, however, easy to see that 
the radicand —duZodyZo can be negative on some surface l^{u,v) (an analysis 
shows that it is OuZq that changes its sign). It means that near this surface 
the linear solution of Eq. becomes inapplicable. Supposing that the GW 
amplitude is everywhere small, 

\Piu)\ « 1 (25) 

and decomposing the radicand at the r.h.s. of Eq. ([24| in series with respect to 
/3, we obtain the equation of this surface: 

I](u,w): {v-u)l3' + 13 + 1^0, (26) 
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which shows that this surface exists at sufficiently large values of the variable 
{v — u) = V2x, i.e., it is far from the boundary: 



\v-u\^ —, 



(27) 



where /3o is the GW amplitude and uj its frequency. 
Near the surfaces v): 



(28) 



whence, putting /3(u) = /3o(l — cos uju), we obtain the equation of the surface 
Yi{u,v), resolved explicitly with respect to retarded time: 



-1 + cr 



Pqoj sin(a;u) + u ' 



(29) 



The ranges in which Eq. p9|) holds are represented by narrow parabolas 
which have almost flat vertices lying near the straight line u — v+a{a = Const). 
Their qualitative behaviour is depicted in Fig. [51 



o 



I! 



u=v+a / 



Figure 2: Ranges of violation of the linear approx- 
imation of Eq. Ii22p 

A detailed image of the vertex of one of the parabolas is shown in Fig. |3l 
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Figure 3: The dependence v(u ) inside one of ranges, 
shown in Fig. [21 

Far from this surface, Eq. ([M)) is integrated, and we obtain the first-order 
correction: 

Zi = - u) e^^C + 2 / C cosh(/3(M)) du (30) 
Jo 

where 

C = / $(w)e'^(") dw. (31) 
Jo 

Since an analytical solution to Eq. (|2T|) cannot be found, while a direct 
application of numerical methods faces considerable difficulties, related to large 
values of the derivatives near singular points, we have used, for a numerical so- 
lution of Eq. (|2ip. the symmetric reflection method, approximating the solution 
behaviour near a singular point with a symmetric parabola. 

Comparing the resulting numerical solution of Eq. (|2ip near singular points 
with the analytical solution ([23| of Eq. ([22]) . we conclude that, despite the 
difficulties indicated above, these solutions almost coincide. An illustration is 
given in Fig. ID 



Figure 4: Comparison of the numerical and analyt- 
ical solutions at small scale inside a range shown in 
\^ the straight line represents an analytical solution, 
the points give a numerical solution. 
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[5] shows a qualitative pattern of the magnetic field strength squared near a 
singular point. As a result, taking into account the nonlinearity of Eq. (|2ip 
near singular2 points reduces, basically, to a cut-off of the lower range (H^ < 0) 
in the graph [T] and to formation of a plateau of the function near zero. Let 
us remark that a similar behaviour of magnetoactive plasma inside a nonlinear 
range was also established in Ref. [1]. 




. 5. MagneticGeld behaviour near singular points 

As a result, taking into account the nonlinearity of Eq. (pij) near singular2 
points reduces, basically, to a cut-off of the lower range {H^ < 0) in the graph 
[T] and to formation of a plateau of the function 77^ near zero. Let us remark 
that a similar behaviour of magnetoactive plasma inside a nonlinear range was 
also established in Ref. |4]. 

Let us also remark that, in a weak gravitational wave (/3 <C 1), the condition 
dvZo w 1 is satisfied, and therefore Eq. ([21]) can be reduced to the form of a 
first-order quasi-linear partial differential equation; 

du(p * (p^^^^ + 4:a^dyip = ^' (32) 

where cp = --duZ. 

It is easy to obtain the formal common solution of this equation, however, 
since it is impossible to find this solution explicitly, it is also impossible to find 
the fmiction Z{u, v) satisfying the required initial and boundary conditions. 
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